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OVERVIEW: This report presents a unified treatment of the density of states for the 
knots and folds of polymer chains. The physical realization of such systems ranges from DNA 
molecules (Taylor) to the microscopic configurations of space-time (Rovelli,Baez). Explicit 
calculations employing this procedure will appear in a subsequent paper. 

The method in brief is as follows: 

1. the starting point is a straight line/chain connecting two boundaries A and B as shown 
in Fig 1 

2. the line is bent at each node, starting at one end and successively moving down the length 
of the chain one node at a time 

3. at each point where the chain crosses itself a " cross-over /under" is assigned 

4. a complete set of configurations, the function space of the knot/fold, is obtained by al- 
lowing all possible folds/crossings (sequences of folds/crossings are characterized by binary 
strings) 

5. a matrix construction is used to characterize the topologically distinct configurations; we 
comment on the implied curvature of the resulting space 
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PART A: folds 



Motivation 

It has been said that the phenomenon of polymer folding is an intractable physics prob- 
lem, yet at first glance the problem appears to be well suited to the machinery of statistical 
mechanics: first choose a model pair potential between joints in the chain, then use the 
total energy in the Boltzmann factor and sum over all possible configurations. Iteration of 
this procedure then follows by suitable adjustments to the assumed pair potential. A diffi- 
culty arises in that the pair potential is likely to be dependent on the entire configuration; 
this is almost certainly due to quantum mechanical influence on the electron orbitals that 
determine the pair potential. Failure of statistical mechanical models suggests the need to 
separate energy and configurational considerations. 

1. Statistical Mechanics 

Models of the self-avoiding walk are frequently used to describe the quantitative proper- 
ties of polymer/protein chain folding (Freed, Pande et al, Amit et al). Typically these tools 
employ standard methods in statistical mechanics. A multi-particle Hamiltonian is set up 



representing a fictitious energy of interaction between two elements in the chain located at 
Xi and Xj . The statistical properties of the chain are then obtained from analysis of the 
partition function 



where /3 acts as an inverse temperature parameter. It is a straight forword procedure to 
evaluate this quantity, but because the interactions must be strong to effectively eliminate 
configurations where the chain crosses itself, standard perturbative techniques are of limited 
value. 

Evaluation of the partition function is not trivial of course, but nevertheless provides a 
complete and systematic solution to the problem. The objection to this approach is that 
the effective/model potential is likely to be unrepresentative. The energy of interaction 
between sites is unavoidably quantum mechanical and thus the interaction will vary not just 
on distance between pairs but also on total chain configuration (Bryngelson 1994). Typically, 
simulations have found the need for what are called non-native interactions (Wallin et al) in 
order to match with experiment. 

To this extent we set energy considerations aside and limit the discussion to configura- 
tional concerns. Exact pair interaction must almost always depend on total configuration, 
hence, non-native interactions. 
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2. The Folding Operator 



Each initial chain is a line on the interval [0, A^]. Each fold describes a bend of ±|, the 
first bend being at the node — 1, the second at N — 2, etc. A binary string [00101. ..Oil] is 
used to charcterize the ± sequence for one complete fold, while a second is used to describe 
the set of crossings as shown in Fig 1. The matrices that perform this operations are as 
follows. 




Figure 1: typical configuration, in this case 3i, the trefoil 



A 



Figure 2: typical bent polymer without crossings 



The definition of the folding operator is that it bends, at right angles, an initially straight 
chain, one site after another. There is one such operator, or matrix in this formulation, for 
each configuration, ignoring self-intersections, and it will be the goal of the analysis to pick 
out those matrices that create allowable configurations. 
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The folding operator is easy to construct and relatively easy to diagonalize but we will 
find limits to the extent of useful information obtained. 
In general the rotation operator takes the form 



cos(^) sin(^) 
— sin(6') cos(6') 



where 9 is the clockwise angle of rotation from the positive x-axis (3 o'clock position). For 

,vr. r 1 
-1 

Define /„ = diag{l, 1, 1) n elements. For a chain of n+1 sites set from (0,0) to (n,0), 



folding at the second to the last site is 



1 



where ri = sir, s = ±1; folding at site third from the end is 



'n-3 



1 

(^2 T2 



etc. 



For example, successive applications of the first four of these matrices yields 
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where a„ = 1 - r„, and TiS2 = SiSsS2r^. 

These are random matrices to be dealt with in the follow-up paper where we discuss the 
distribution of eigenvalues. Preliminary numerical studies suggest allowed structures tend to 
cluster in groups. Such behavior has been reported in previous studies (Balafras and Dewey, 
Moret et al) and has possible application to the so-called Levinthal Paradox (Karplus, Dill 
and Chan). 

The Levinthal Paradox says that a large polymer would take eons rather than millisec- 
onds to fold if it did so by randomly sampling accessible states. What this suggests is that 
the folding process (and knotting as well) is accomplished through some globally determined 
energy potential and is not a locally driven phenomenon (Wallin et al). 
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Part B: Knots 



Motivation 

The difference between a Icnot and a fold in this formulation of the problem is the nature 
of the sequence of crossings: a particular sequence of crossings may, or may not, actually tie 
the chain, such that it becomes possible, or not, to pull the two ends arbitrarily far apart, 
returning the chain to the initial straight line. 

How to do this is dealt in part by the following scheme. Each knot is thought to be de- 
scribed by two dimensions, position (or crossing) and time (or path length). The distinction 
as suggested in Fig 3 is made by first labeling all crossings successively until each site has 
a number, followed by assigning to each site a "time" value as a full circuit of the knot is 
traversed. 

[Aside: For the first few knots all crossings are labeled before the return path is be- 
gun; these are called sequential. However, beginning with Ty a previously labeled site is 
encountered before all seven sites have been labeled; these are called nonsequential. This is 
something like Euler's Konigsberg Bridge problem. This distinction is made clear in Figs 4 
and 5 where we show how to reduce 63 to a simpler circuit, thus making it clear that the 
knot is sequential, and also show that 7j is necessarily nonsequential.] 



1. Knot Matrix Designation 

For a knot with N crossings: 

a) each crossing is an element in a sequence of positions 

b) the n^^ position element is described by the matrix 



Tn{t) 



n-1 



T{t) 



'Af-n-1 



where 



T{t) 



e 

1 



iat 



for cross-over and for cross-under; a to be determined. 

c) the time sequence for a given knot is the product [T] = TaTbTcTd... where abed... are 
successive crossings visited on a complete path; for example for 61 the [T] product is 
T5TQT1T2T3T4TQT5T4T3T2T1 abbreviated to [561234.654321]. This procedure is not unlike 
the group theory method of Artin/Burau (Kauffman, p. 86) except that the operators are 
2d and label crossings instead of braids (Birman and Brendle). Also it appears that these 
matrices do not obey the Yang-Baxter relations (Baez and Muniain). 

d) as in Fig 3 each operator is 2d, i.e. Ta{n) where (a,n)=(space,time)= (position, path length). 

The matrices r, act at each crossing something like raising and lowering operators in 
field theory, and create a "time evolution" of the initial, ground state (the identity matrix), 
carrying it through a sequence of intermediate states. The resulting matrix is a measure of 
the curvature of the knot to the extent that the configuration is not returned to the identity 
matrix for the closed contour. 
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Figure 3: geometry of the matrix Ti{n) 




Figure 4: construction showing that 63 is sequential 



2. Topological Subspaces; Curvature 

In order to construct the representative matrices for each knot we need to assign a value 
to the time variable and this is accomplished by requiring that configurations isomorphic 
to the unknot give the identity as shown in Fig 6. For numerical purposes we can choose an 
arbitrary real number q and require that the second pass through the crossings act in reverse 
time; for example for 4i the time factors are g, q"^, q^, q'^, g^, q^, g^, q. 

The resulting matrices act as a representation of the curvature of the knot in the sense 
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Figure 6: the identity operation for N=3 



that a complete path around the knot is analogous to the closed curve in the computation 
of parallel displacement of a vector defined on a smooth manifold. The resulting matrix 
bares some relation to the curvature tensor in differential geometry; in quantum field theory, 
this quantity is a relation of the Wilson Loop, without the extra machinery for field theory 
purposes. 

The table below gives the knot, the path and the resulting matrices for the figures shown 
in the Appendix. 
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[Aside: I gave up trying to get LATEX to put the figures where I wanted. The Appendix 
included here contains the relevant knots of the table; the drawings, inexpert as they are, 
indicate the sequence of crossings used to obtain the matrices below.] 
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Part C. Summary; Follow-up 

The methods presented above provide a means for constructing folds and characterizing 
disjoint topological spaces defined by each distinct knot. In a following paper we will fur- 
ther discuss the properties of these subspaces and how they might be used in the partition 
function for statistical purposes. 

The broader implications of this work were not apparent at first. Initially, I had thought 
of the [T] product in the context of polynomial knot invariants, guided by an analogy with the 
Artin/Burau group representation. Only afterward did it appear that the [T] product was 
not unlike the Wilson Loop which is used in quantum field theory: this expression could be 
used to derive the Jones polynomial (Witten, Kauffman, Baez and Muniain). However, the 
matrices appear to contain more information than the polynomials; for example, subspaces 
apparent in the matrix suggest symmetries among subsets of crossings and might act to 
better characterize knot classes. Furthermore, close relation, at least in appearance, to the 
partition function suggests additional associations. 

And yet, the crossings are, after all, not physical. The internal symmetries of knots 
have long been a matter of considerable interest (Hoste et al, Gruenbaum and Shepard) 
and some of the information pertaining to these symmetries is apparent in the mixing (or 
non-mixing) of crossings as indicated in the knot matrices. For example, one would expect 
knots 3i, 5i, 7i, ... to be diagonal, as they are (see below). 

The degree of mixing of the crossings, as basis vectors, is in some sense a measure of the 
complexity of the knot. The entropy involved in such configurations has recently received 
considerable attention (Baiesi et al). 

Some open questions to be considered in the follow-up paper: 

1. it is not clear how to relate the matrices to invariant polynomials, such as Alexan- 
der/ Jones/Kauffman, nor how, if at all, to construct skein relations for the [T] products 

2. one change in methodology, possibly, periodic boundary conditions suggest that the r 
matrix should, for the last crossing on the path be 

" ... e*"* ' 
1 ... _ 

instead of 

/iv_i ... 

e*"* 

1 

This change in procedure would assure the diagonal behavior of the simple odd numbered 
knots. 

3. the notion of curvature of a knot is admittedly vague, unlike that of, say a 3-d manifold 
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Figure 9: 6i 62 knots 




Figure 10: 63 knot 
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Figure 12: 74 
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Figure 14: 
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